An electron in the presence of multiple zero range potentials and an external laser 
field — exact solutions for photoionization and stimulated bremsstrahlung 
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The method of zero range potential (ZRP) for one-electron problems is reviewed. In the absence 
of an external electromagnetic field, the notion of a ZRP is introduced from different points of view 
and for an arbitrary dimension of space. Then, three-dimensional problems of motion of an electron 
in the field of several ZRPs and laser radiation are studied. Exact wave functions for the processes 
of photoionization and stimulated bremsstrahlung in the presence of the laser field of an arbitrary 
pulse shape are obtained in the form of one-dimensional integral representations. 
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I. INTRODUCTION 

The zero range potential (ZRP) method is widely 
employed by physicists, mathematicians, and chemists 
P) S, S 0, Si because it usually allows one to obtain ex- 
act solutions^ for a broad variety of model systems. The 
term ZRP may sound bizarre for beginners because a 
potential is usually represented by a non- vanishing func- 
tion in some region. An example of the "function" be- 
ing nonzero at a single point ro is the Dirac ^-function, 
(5(r — ro). Thus, the ZRP has various other names such 
as, e.g., "point interaction," "delta interaction," "contact 
interaction." But the (5-function is not the only example. 
In general, an arbitrary "function" concentrated at rg 
can be represented as ^„ c„5^"^ (r — Tq) [S]. However, all 
these objects are not ordinary functions; they are distri- 
butions and require special mathematical care [6|. (We 
will return to this way of tackling a ZRP in Sec. IIIDI ) 

The term ZRP should not be taken too literally. In 
fact, the ZRP method should be understood rather as a 
schematic or axiomatic description of the interaction of 
a particle with an object, which we call the ZRP. Note 
that there is no standardized definition of the ZRP; the 
definition varies for different problems. The essence of 
the ZRP method is to define what we mean exactly by 
the ZRP in a given situation. In the current paper, we 
pursue this point of view on the ZRP and demonstrate 
different approaches to define such a notion and inter- 
relations among these definitions. Such a presentation 
should be beneficial especially for neophytes. 

Kronig and Penney ^ were the first to employ ZRPs. 
They formulated the well known Kroning-Penney model 
of solids. A few years later, Bethe and Peierls [8j and 
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Thomas [9| used a single ZRP in the theory of a deuteron. 
The first application of ZRPs in atomic physics was made 
by Fermi [IQ] . He also considered neutron scattering in 
substances containing hydrogen by means of ZRPs 
The first rigorous mathematical treatment of the notion 
of a ZRP was given by Faddeev and Berezin [12]. Re- 
garding further development of the method of ZRP see, 
e.g., Refs. [1, H, [1] and references therein. 

The method of the ZRP entered strong field physics 
when analytical results on ionization of a weakly- 
bounded electron (negative ions), i.e., bounded in the 
field of a ZRP, were reported [3 H El- The exact 
three-dimensional solution to the problem of ionization 
of a particle in a ZRP by circularly polarized laser ra- 
diation was found in Refs. Afterwards, ioniza- 
tion in a elliptically polarized monochromatic laser field 
has been solved exactly in Refs. [H, [l^. Stimulated 
bremsstrahlung radiation in the presence of a ZRP was 
studied in Refs. [itIjUO]- The exact solution for the one- 
dimensional Kronig-Penney model in a monochromatic 
laser field was achieved in Refs. [2l|. High-harmonic 
generation from a model atom with a ZRP was inves- 
tigated in Ref. [l^l- The ZRP method was used to 
model quasimolecular ions in a circularly polarized laser 
field in Ref. [1^. The exact solution for the prob- 
lem of laser-assisted scattering from a one-dimensional 
(5- function potential was reached in Refs. [13, [131 • De- 
cay of a weakly bound level in a monochromatic elec- 
tromagnetic field and a static magnetic field was stud- 
ied in Ref. [1^. Photoionization of a negative ion in 
a bichromatic laser field was analyzed in Refs. [H, [13| ■ 
Recently, the ZRP method and its generalizations have 
been remarkably useful in understanding rescattering, 
above threshold ionization, and hig h harmonic genera- 
tion m, [H, [H [M [13, [11, [13, [If. A more complete 
list of applications of ZRP in strong field physics can be 
found in reviews [H, [s^l • 

Despite intense research in this direction, the problem 
of photoionization of a weakly bounded electron by the 
laser field of an arbitrary vector potential, A(t), has not 



2 



been solved exactly. Not only do we solve such a problem 
in the current paper, but we also obtain solutions for pho- 
toionization as well as stimulated bremsstrahlung of an 
electron bounded in the field of several ZRPs. The latter 
one is not only of an academic importance. Processes in- 
duced in molecular systems by an external laser field is a 
hot topic in strong field physics (for review see, e.g., Ref. 
[3^). Evidently, Coulomb forces ought to play an impor- 
tant role in these processes. Nonetheless, identification 
of features caused by the Coulomb interactions is a chal- 
lenging task in general case. Our solution would allow 
to answer a complementary question: "which processes 
are not due to the long range part of the Coulomb forces 
?" Indeed, the system of an electron moving in the field 
of many ZRPs is perhaps the simplest molecular model 
without long range forces. 

It is noteworthy to mention that the ZRP method is 
not only fruitful for modeling a broad variety of so-called 
weakly bounded systems, but also it has been utilized to 
solve approximatelythe Schrodinger equation for some 
class of potentials |39l. [40|. Furthermore, the 5- function 
as a potential term in the Schrodinger equation has the 
following interesting property: the Hamiltonian in this 
equation has the minimum lowest eigenvalue among all 
potentials of a given "area" ,41]. Hence, the convergence 
of variational calculations may be improved by includ- 
ing solutions of the corresponding ZRP problem into the 
variational basis set (see, e.g., Ref. Q). 

The structure of the current paper is as follows. As it 
was mentioned at the beginning of the Introduction, the 
notion of a ZRP is multifaceted; an attempt to manifest 
its richness by presenting different definitions of a ZRP 
is made in Sec. HT] Presenting the matter of Sec. [Tll we 
often sacrifice mathematical rigor and technical details in 
favour of a clear and physically intuitive presentation; ad- 
ditionally, we present references to literature where more 
detailed discussions can be found. In general. Sec. |TT] 
can be regarded as an overture to Sec. IIIIl which is the 
main part of this paper. Problems of photoionization 
and stimulated bremsstrahlung of a one electron system 
in the field of ZRPs are discussed and solved exactly in 
Sees, nil Al and IIIIBl correspondingly. 

Finally, for the sake of clarity, we mention a few words 
regarding notations and conventions used thorough the 
paper. Let d denote a dimension of the space. Many 
equations are valid for an arbitrary dimension of the 
space; thus unless stated otherwise, the positive inte- 
ger d is not specified. Bold letters (e.g., r, k, p) de- 
note d-dimensional vectors; r = |r| - the absolute value 
of (i-dimensional vector r. S{r) is the d-dimensional 
Dirac (5-function. N denote the number of ZRPs, and 
unless stated otherwise its value is not fixed. Vectors 
Rj (j — 1,2,..., N) specify the locations of the ZRPs. 
Atomic units, h = me = |e| = 1, are used throughout. 



II. A ZRP WITHOUT THE PRESENCE OF A 
LASER FIELD 

A. Physically intuitive method 

Probably most straightforward and intuitive way of in- 
troducing a ZRP is to replace the Schrodinger equation 
inside the ZRP by a certain boundary condition on the 
wave function at the center of the ZRP. This boundary 
condition can be postulating after analyzing solutions of 
the Schrodinger equation for a "shrinking" potential. Let 
us present an example of such a description \^ [4^ . 

Assuming the three-dimensional space, we consider a 
particle in an s-state in the presence of a spherically sym- 
metric potential well defined as rj < b, where Vj = r — Hj 
is the distance between the position of the particle r and 
the center of the well Rj and b is the radius of the well. 
The wave function of the particle outside the potential 
well has the following form 

V'(r) = Cjexp{-ajrj)/rj, (r^ > b) (1) 

where Cj is the normalization constant and aj = \J'2.1p^ 
Ip is the binding energy (ionization potential). Now we 
contract the width of the well {b 0) and simultaneously 
increase its depth {Uq oc 1/&^) such that Ip as well as aj 
remains constant. Hence, the potential well is replaced 
by a boundary condition at Rj-: 

V'(r)'^^°c,[l/r,-a,+0(r,)]. (2) 

Having done this limiting procedure, we can define the 
ZRP as follows: a wave function tp{r) of a particle moving 
in the field of a ZRP located at Rj is a solution of the sta- 
tionary Schrodinger equation, — AV'(r)/2 = Eip{r), sub- 
jected to boundary condition 12]). Note that parameter 
aj does not depend on the energy of the particle E; aj 
is to be regarded as merely a characteristic of the ZRP 
and 2t: /aj is sometimes called the depth of the ZRP Q 
or the renormalized coupling constant (see Sec. Ill Bp . 

One might have a methodological concern regarding 
boundary condition ([2]): can a wave function have a sin- 
gularity in quantum mechanics? Indeed, the postulates 
of quantum mechanics do not forbid for a wave function 
to be singular as long as the probability of finding the 
particle in a volume V, 

Pr(r e = / |i/.(r)pd3r, (3) 
Jv 

which is a measurable quantity, is continuous. Since the 
integral in Eq. ([3]) exists if y is a small ball centered 
at Rj, wave functions satisfying boundary condition ([2|) 
are allowed. Being more rigorous, we recall the first pos- 
tulate of quantum mechanics: that wave functions are 
described by normalized vectors of a Hilbert space. The 
Hilbert space under discussion is ^2(1^'^) - the set of func- 
tions 0(r) such that J \(j){r)\'^d'^r is finite; if 1/ = R'^, then 
Eq. ^ is the definition of the norm of such a Hilbert 
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space. In other words, ^2(R ) is a set of bound states. 
Thus, it is noteworthy to recall the fact that scattering 
states, which are solutions of the stationary Schrodinger 
equation subjected to outgoing or ingoing wave bound- 
ary conditions at infinity, are not elements of the Hilbert 
space ^2(1^'^) simply because they are not square inte- 
grable. Mathematically speaking, the scattering states 
are distributional eigenfunctions of the Hamiltonian (re- 
call that they are "normalized" to the 5- function). 

The presented limit transition can be performed in the 
one- and two-dimensional cases. The boundary condition 
for a wave function ^'{x) in the one-dimensional case is 
given by 



dx 



x=Xj+0 



dx 



(4) 



where Xj denotes the position of the ZRP; the boundary 
condition for the two-dimensional case reads 

V'(r) C.J [In Tj + ln(a, /2) + C + 0{r] In )] , (5) 

where C = 0.5772 . . . being the Euler-Mascheroni con- 
stant. 

Further important generalizations arc possible. 
Boundary conditions ([3]), and ([5]) remain unaltered 
even if we add a potential without a singularity at Rj 
into the stationary Schrodinger equation. In the case of 
many [N) ZRPs, we have N boundary conditions for 
J = 1, 2, . . . , iV (constants Cj as well as Uj may differ for 
different j). 

The stationary two-center problem within the ZRP 
method [N = 2) has been pioneered by Smirnov and 
Firsov [4^; afterwards, the general solution for the prob- 
lem of a particle motion in a combine fields of several 
ZRPs has been found [l^H^. Moreover, there have been 
obtained a broad variety of Green's functions for such a 
problem (see, e.g., Refs. [U IH, H^). Regarding cur- 
rent status of the usage of the method of ZRP to model 
molecular systems see, e.g., Refs. [50]. 

Different boundary conditions also have been obtained 
and applied to more elaborate cases, e.g., such as a par- 
ticle in the presence of an external static magnetic field 
psi . [sij , the time-dependant formulation of detachment 
of an electron from a negative ion in collision with an 
atom [5I], and many others. A more complete list of 
applications can be found in Refs [3, [1] . 



B. Self-adjoint extensions of a Hermitian 
Hamiltonian 

This section is an effort to build the bridge between two 
stylistically different points of view on the same object - 
mathematicians' and physicists' way of defying a ZRP. 

According to the postulates of quantum mechanics, 
measurable quantities are represented by a self-adjoint 
operators, which act on a Hilbert space of wave functions 
of a given system. Note that there is a subtle difference 



between a Hermitian and self-adjoint operator. Usually, 
this difference is not important for physicists; however, 
it plays a crucial role in the theoretical formulation of a 
ZRP. 

Let us briefly clarify the difference. A linear opera- 
tor A, acting on a Hilbert space J^, is a linear mapping 
A : D{A) J^, where D{A) is a domain of the operator 
- a set of elements on which the action of operator is 
defined. (In general, the domain D{A) does not coincide 
with the entire Hilbert space Jif.) The adjoint opera- 
tor At of an operator A, whose domains are D{A) and 
D{A^) correspondingly, is such that the following equal- 
ity is satisfied 



(it/,5) = (/,i5), V/,5ei?(i). 



(6) 



The operator is sclf-adjoint if A'' = A, i.e., the action of 
coincides with the action of ^, and D{A) C D{A''). A 
self- adjoint operator A is a Hermitian operator if D{A) — 
D(At). Note that there is no distinction between self- 
adjoint and Hermitian operators in a finite dimensional 
Hilbert space J^f because operators are given by matrices 
and the operation of multiplication of a finite matrix by 
a vector is defined for all vectors from , viz., D{A) = 
D{A'') = ^ for an arbitrary matrix/operator A. 

Now we summarize two important results from the the- 
ory of self-adjoint operators that are widely employed in 
physics often without being referred to. The purpose 
of such a summary is to clarify the postulate of quan- 
tum mechanics, i.e., to accentuate the importance for the 
Hamiltonian of a quantum system as well as other ob- 
scrvables to be represented by self-adjoint operators and 
not by Hermitian ones. The first result is the Hilbert- 
Schniidt theorem [ssj : if is a compact self-adjoint op- 
erator, then its eigenfunctions form an orthonormal basis 
for the Hilbert space. Second, Stone's theorem [s^, H^: 
if U{t) is a strongly continuous unitary group (which de- 
scribes the evolution of a quantum system), then there is 
a unique self-adjoint operator H (the Hamiltonian of the 
given quantum system) such that 

^§-Uit)f = HU{t)f, V/ e D{H), (7) 

this equation is formally expressed as U(t) = exp(— ii/t). 
Physically speaking, if the time-independent Hamilto- 
nian H is a self-adjoint operator, then Stone's theorem 
guarantees the unique solution of the time-dependent 
Schrodinger equation such that the normalization of the 
wave function is constant in time. (Regarding the gen- 
eralization of Stone's theorem to the case of a time- 
dependent Hamiltonian see, e.g.. Sec. X.12 of Ref. 56].) 

An important fact is that a Hermitian operator can 
be transformed into a self-adjoint one by redefining the 
domain of the operator, e.g., by means of imposing ap- 
propriate boundary conditions on the wave functions on 
which the operator acts. This procedure is an object of 
study of the theory of self-adjoint extensions of a Her- 
mitian operator. We shall not discus this theory in the 
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current paper. There are tutorials oriented for physicists 
that not only explain this theory, but also illustrate it on 
many physically interesting examples [HI, [s^, [H, [sl, [sO] ■ 
Regarding the review of the theory of extensions and 
exactly solvable models see, e.g., Ref. [gj]. Rigorous 
discussions of the theory can be found in mathematical 
textbooks [13, H^- (There are also textbooks oriented 
for applications of functional analysis [53, 56].) 

We describe the application of the theory of self-adjoint 
extensions of a Hermitian Hamiltonian to the problem of 
ZRP in a nutshell (for further details and more rigor- 
ous formulations see Ref. Q as well as tutorial [6C|). 
This treatment of a ZRP by means of the theory of self- 
adi oint extensions was pioneered by Faddeev and Berezin 
Assuming that a single ZRP is located at the ori- 
gin, we denote the Hamiltonian for a particle in the 
field of the ZRP by Hzrp and d = 1,2,3 - the di- 
mension of the space. A priori fact is that the Hamil- 
tonian Hzrp must coincide with the free particle Hamil- 
tonian, Hq = —A/2, in the space without the origin, 
M.'^ \ 0. Therefore, let us restrict the domain of Ho to 
'^o°°(I^'' \ 0) - a set of infinitely differentiablc functions 
that vanish outside a compact set which does not con- 
tain the origin (e.g., let d = 1 then '^q°°(R\ 0) consists 
of infinitely differentiable functions fabcd{x) that vanish 
for X € (— oo, —a) U {—b, c) U (d, -|-oo), where a, 6, c, and 
d are arbitrary positive numbers). One can verify that 
Hq defined on such a domain is a Hermitian but not self- 
adjoint operator. Applying the theory of extensions, one 
concludes that for d — 2,3 there is a one-parameter fam- 
ily of self-adjoint extensions of the Hermitian operator 
Hq indexed by a "renormalized coupling constant" [and 
defined on an appropriate subset of J!f2{M.'^)-] These self- 
adjoint operators are to be postulated as definitions of 
Hzrp in the two- and three-dimensional spaces. How- 
ever, the one-dimensional case is special because the Her- 
mitian operator Hq has a four-parameter family of self- 
adjoint extensions, which physically means that there are 
not only (5-interactions [this alternative name for a ZRP 
comes from the fact that heuristically speaking, a one- 
dimensional ZRP can be represent as the term —^S{x) 
in the Hamiltonian, Hzrp — —A/2 — fj,6{x)], but also 
new types of point interactions such as (^'-interaction 
[which has a heuristic form —fi6'{x)] and powers of the 
5- interaction [gl, [6^ . Note that these new interactions 
cannot be introduced by means of the physically intuitive 
method discussed above. 

To draw the connection between this section and sec- 
tioning we mention that the operator Hq defined on the 
set of functions that satisfy boundary conditions ([1]), ([1]), 
and ([5]) is indeed self-adjoint. Unfortunately, book ;3i| as 
well as other physicists' publications demonstrates a lack 
of understanding of the last fact - Hq is demonstrated 
to be a Hermitian operator (see page 6 of Ref. (H), but 
no comment is made regarding its self-adjointness. This 
is not remarkable since physicists rarely distinguish be- 
tween Hermitian and self-adjoint operators. Needless to 
say, such a distinction is vital as far as the notion of a 



ZRP is concerned. 

Now we clarify why Hamiltonians that contain the S- 
function as a potential term [e.g., Hzrp = —A/2— fiS{x)] 
are not legitimate, and thus they are to be regarded 
as heuristic expressions. The quantity fJ-S{r), where 
/i G M, is not an operator acting on the Hilbert space 
^2{^ ) since for an arbitrary /(r) e if2(K''), the inte- 
gral / \ S{r)f{r)\'^d'^r does not make any sense. From this 
point of view, we want to emphasize that the approach 
of self-adjoint extensions of the free particle Hamiltonian 
does not assign any potential to the interaction of a par- 
ticle with the ZRP, viz., the method of self-adjoint ex- 
tensions gives us the "potential-free" interpretation of a 
ZRP. 



C. Renormalization Method 

A variety of renormalization techniques, some being 
similar to renormalization in quantum field theory, have 
been used in the problem of ZRPs [H, [H, 113, US', 69| (see 
also Ch. X.ll of Ref. (HI). In this section, we advocate a 
renormalization method, which combines advantages of 
known renormalization techniques. This method will be 
adapted to the time-dependent case in Sec. IHII 

We introduce our renormalization method on the ex- 
ample of N ZRPs that are located at {Rjjj^^, and the 
dimension of the space d is not specified. We start from 
the time-independent Schrodinger equation 



(8) 



where V is the heuristic potential, which represents the 
ZRPs, and whose coordinate representation reads 

N 

{r,\V\r2) ^-S{r,-r2)Y.fi,S{r,-R,). (9) 

Projecting Eq. ^ onto the conjugate of a momentum 
eigenstate (k|, and using the equality (k| Rj) — exp(— ik- 
Rj)/(27r)''/^, we obtain the Schrodinger equation in the 
momentum representation 

N 

{E - kV2) (k| f^J exp(-zk • R,), (10) 



where 



P,=f,,{R,\i,)/i2nf/' 



(11) 



being a constant. Solution of Eq. PH]) depends on the 
sign of the energy E: ii E < 0, then there is a bound 
state solution 



N 



(k|V)=E 



l3j exp(— «k • Rj) 
ky2+\E\ ' 



(12) 
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which will be an element of ^2(1^'') after Fourier trans- 
formation. If E is positive, then taking into account that 
(p^ — k^)(5(p— k) — 0, we obtain scattering state solutions 



(k| ^) - S{p - k) 



N 

E 



(3j exp(— ik • Hj) 
k^/2-E±iO 



(13) 



where = E and ±iO correspond to outgoing and 

ingoing wave boundary conditions, respectively. Equa- 
tion (fT3|) clearly illustrates the fact mentioned in Sec. 
Ill Al that scattering state solutions are not elements of 
^2(R'') ^ they are distributions. 



Now we scrutinize bound state solution (|T2|) . Flourier 
transforming wave function p2p . we obtain the wave 
function, ip{r) = (r| ■0), in the coordinate representation 
in the following compact form (for an arbitrary spacial 
dimension d) 



N 



{d-2)/2 ' 



(14) 



where Bj are constants proportional to (3j, k — ^y2\E\ = 
y/2Ip, and K^{x) is the modified Bessel function of the 
second kind (or the Macdonald function [tO^). We point 
out that the ZRP in higher dimensions {d > 3) is also of 
interest (see, e.g., Refs. [13, IZll)- Let us itemize inter- 
esting special cases of Eq. (fT4| 



N 



Bj exp(— kIxj I), if d = 1, 



N 



^BjKoinrj), if d = 2, 



(15) 



3 = i 



N 



Bj exp(— Krj )/rj, if d = 3. 



In the case of a single ZRP {N — 1), one can readily 
verify that wave function ()14p indeed satisfies boundary 
conditions ([4]), and ^ for d = 1,2, and 3, corre- 
spondingly. In the case of many ZRPs, following Refs. 
psi . [46j (see also Ch. 3 of Ref. we impose the 

boundary conditions on wave function (fT4| at the points 
{r = Ilj}jLi and obtain a homogeneous system of A^ 
algebraic equations, from which the coefficients {Bj}jL^ 
and K can be obtained. 

Now we shall return to the definition of the parame- 
ters l3j [Eq. (dJ)]. Since (Rj|V') = V^R-j) and i/'(Rj) 
are singular for d ^ 2, we conclude that the parameters 
(3j contain divergences in the case of c? ^ 2. However, 
we have not used the freedom of choosing values of the 
coupling constants /ij. Hence, we interpret the values of 
ipCRj) as some infinite numbers, and thus set fj,j to be 
equal to corresponding infinitesimals such that (3j G M. 
Summarizing, we conclude that the normalization con- 
stant of the wave function {Bj cx (3j) is renormalized 



within such a method. Note that the one-dimensional 
case is special because renormalization is not needed. 

The renormalization technique is generalized to de- 
scribe a bound state with nonzero angular momentum 
in Appendix [Al 



D. Miscellaneous methods 

In this section, we list other methods used to define 
the notion of a ZRP. Nevertheless, this list will not be 
exhaustive, and we refer the reader to monograph Q for 
other approaches. 

In Sec. Ill C[ we constructed the wave function of a 
particle in the field of ZRPs by employing heuristic po- 
tential ([9]) with infinitesimally small coupling constant fi. 
However, one may still hope that it is possible to repre- 
sent consistently a ZRP by means of the (5-function with 
a non- vanishing coupling constant and additionally avoid 
the problem that the (5-function does not define an op- 
erator on the Hilbert space (see Sec. Ill Bp . As we have 
seen (e.g.. Sec. Ill A[) that the wave function or/and its 
derivatives has a singularity at the point where a ZRP is 
located. Since the classical calculus is not applicable to 
the wave function at the singular point, we shall treat it 
as a distribution to have the properly defined operation of 
differentiation. If we follow such a path we immediately 
run into another problem because the multiplication of 
distributions is not uniquely defined, and the stationary 
Schrodingcr equation contains the product of the wave 
function and the potential, which are distributions. To 
illustrate the troublesomeness of multiplication of distri- 
butions, we recall the Schwartz counterexample, 

= 0-l/x= {d{x) ■x)-l/x ^ 6{x) -{x-l/x):^ d{x), 

i.e., there is no associative and commutative operation of 
multiplication consistent with multiplication by multipli- 
ers. The problem of multiplications of distributions is an 
active area of research in modern mathematical physics; 
some approaches to this problem can be found in Refs. 
[64, 65, 72j, [tI, [tI]. Summarizing, we say that for the 
construction of quantum mechanics with singular poten- 
tials, it might not be sufficient to regard functions cor- 
responding to singular potentials as distributions due to 
onerousness of constructing a sufficiently large associa- 
tive algebra of distributions. 

Nevertheless, there have been successful attempts to 
construct axiomatically suitable algebras - Shirokov's al- 
gebras [zi. III, [77I [tsI. [79I, and they have been used to 
tackle the ZRP [lalsin^- Antonevich [Hi has also 
employed Colombeau's algebra of new generalized func- 
tions [731 to the problem of ZRP. As a historical note, we 
mention that in the fifties, Bogolubov suggested that a 
necessary attribute of any future local theory would be 
a modification of the subtractional formalism associated 
with the multiplication of distributions. Furthermore, 
there is a strong connection between regularization in 
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quantum field theory and the problem of multiplication 
of distributions 

[Mil- 

In Ref. [s^, a scalar product of some distributions 
was introduced and used to define the one-dimensional 
ZRP as a singular potential bypassing the problem of 
multiplications of distributions. In some sense, this work 
is in spirit of Shirokov's algebras. 

The renormalization technique (Sec. Ill Cp may also 
suggest that an elegant and satisfactory picture of ZRP 
could be achieved if the notions of "infinite and in- 
finitesimal quantities" were formalized. This point of 
view leads to the nonstandard analysis, which can be 
looked upon as a rigorous mathematical realization of 
Newton's and Leibniz's ideas of infinitesimal quantities 
without running into contradictions. Non-standard anal- 
ysis was introduced in the early sixties by Robinson 
[87| . (Non-mathematicians may find a comprehensible 
introduction to the nonstandard analysis in book [s^].) 
In Refs. [S^, [sO], the ZRP was quite organically de- 
fined within nonstandard analysis. Moreover, there have 
been attempts to formulate quantum mechanics in a non- 
standard Hilbert space 

Fermi [10] indirectly proposed to replace the expression 
V = — /i(5(r) in the three-dimensional case by 



27r _ , 9 2tt ^ , _ 
Vf = —S{r)—r = —6ir)V ■ r, 
K or K 



(16) 



where as previously n ~ 



and we assumed that 



there is a single ZRP located at the origin. New heuris- 
tic potential (|16p has been called the Fermi pseudopoten- 
tial, and it has been widely exploited in physics (see, e.g., 
Refs. [13, [63, [11] and references therein). A main ad- 
vantage of the Fermi pseudopotential compare to heuris- 
tic potential ^ is the fact that the coupling constant, 
2'k/k, is a real non-vanishing number. Nonetheless, the 
pseudopotential is still not an operator on the Hilbert 
space, similarly to V (see Sec. IIIB|) . Additionally, it el- 
egantly avoids the problem of multiplication of distribu- 
tions: since Vp acts on a wave function with the asymp- 
totic behavior ([2]) near the origin, we obtain 

Vf [l/r + K + 0{r)] = 27r(5(r). 

However, the situation is opposite once phenomena in an 
external electromagnetic field are considered: Let H = 
— V^/2 -I- Vf denote the Hamiltonian (in the coordinate 
representation) of a particle moving in the field of a ZRP 
modeled by the pseudopotential ^B)) . As it was pointed 
out in Ref. [9^, using the minimal coupling procedure, 
— iV — > — iV -I- A, to "turn on" the electromagnetic field 
defined by the vector potential A — A(r, i), we update 
the form of the Hamiltonian H to 

H^-[-iV + Af + —S{r)[^iV + A]-r. (17) 
2 K 

If we assume that the Hamiltonian (fT7|) acts on wave 
functions that obey boundary condition ^ (or any sin- 
gular boundary condition at the origin), then the last 



term of Eq. ((TT)) indeed represents the multiplication of 
two distributions - the J-function and the wave function. 
Note that the standard symbolic relationship, ri5(r) = 0, 
is valid only when the functional rS{r) acts on an in- 
finitely differentiable function. (The latter fact neverthe- 
less was not recognized in Ref. [OJ]; besides, the existence 
of the last term of Eq. (fTTl) was not even mentioned in 
Refs. 16, 17, 18, 19, 20, 36, 37].) Therefore, we do not 
use the Fermi pseudopotential in the current paper; in 
fact, we firmly discourage one from utilizing it because, 
as it was seen in the previous illustration, the Fermi pseu- 
dopotential is "too" ad-hoc to be useful for the problems 
under scrutiny in Sec. IIIII 



III. 



ZRPS IN THE PRESENCE OF A LASER 
FIELD 

A. Photoionization 



Similarly to scattering states (see Sec. Ill Al) . a wave 
function that describes photoionization is not normaliz- 
able, i.e., it is not an element of the Hilbert space ^2{^'^)- 
Hence, we cannot use the theory of self-adjoint extensions 
(Sec. Ill Bp to obtain boundary conditions, which should 
be satisfied by the wave function. Furthermore, the way 
of defying the ZRP as a limit of some continuous poten- 
tials is cumbersome analytically; nevertheless, one could 
attempt to use numerical solutions of the Schrodinger 
equation with a time-dependent Hamiltonian to find such 
a definition. However, there is an easier way: A straight- 
forward generalization of the renormalization technique 
presented in Sec. Ill CI allows us to define the problem of 
photoionization of a weakly bounded system (bounded in 
the field of ZRPs). 

From now onwards, we solve the time-dependent 
Schodinger equation 



^ + V + p-A{t) + ^A\t) 



l*W>,(18) 



where V is the heuristic potential, which represents the 
ZRPs, given by Eq. ([5]); A{t) is an arbitrary vector po- 
tential. A crucial step to solve such a problem is to make 
the following substitution into the Schodinger equation 



\^it))=eM-^^t)mt)) 



(19) 



Let the substitution (|19p be called the Floquet ansdtze; 
the reason of such a name will be discussed later in this 
section. The Schodinger equation for the wave function 
]$) reads 

^^^mt))-H,Xt)Mt)), (20) 
H,{t) = pV2 + V + P- A{t) + A^t)/2 - e, 
Partitioning the total Hamiltonian 

HS)^Hv{t) + V, (21) 
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we can write down the Lippmann-Schwinger equation in 
the post form for the propagator (see, e.g., Refs. [4, 

U,{t,to) = Uv{t,to)-i I Uv{t,t')VUe{t',to)dt', (22) 



where the propagators and Uy correspond to the 
Hamihonian iJ^ and Hy respectively, viz.. 



He^v{T)dT 



(23) 



where T being the time ordering operator. The exact 
expression of the propagator Uy is known 



Lfvit, t') = exp[ze(t - t')]UvG{t, t'), 



(24) 



where UvG{t,t') being the velocity gauge Volkov propa- 
gator 

KvGir,t\v',t') = {v\UvG{t,t')\v') = 

0{t-t') \ i 

exp 



A{T)dT 



(25) 



Since the problem of ionization is being considered, we 
introduce the following initial condition: the wave func- 
tion |$(t = — oo)) = lip) is a bound state [see, e.g., Eq. 
(fT4|l ] . Substituting the following equalities 



mt))^U,{t,~^)\ij), hm (r\Uvit,to)\^P)=0, 



into Eq. (|22p . we obtain the equation for the wave func- 
tion, $(r,t) = (r|$(i)). 



N 



^r,t) = 5^G,(r), 



(26) 



G,(r) 



oie(t-t') 



;KvGir,t\Il,,t')f,it')dt'i27) 



where fj {t) is a time-dependent analogue of Pj [Eq. pTjl ] 

/,W=*M,<i>(R„t)/(2^)'^/2. (28) 

In the rest of the current section, we focus on the 
three-dimensional case. Cases of other dimensions can 
be tackled similarly. The asymptotic expansion of Gj (r) 
for small rj is found in Appendix [B] [Eq. (|B4p ] . Now we 
are in position to postulate boundary conditions in the 
case of the presence of a laser field. We shall be guided 
by the two following principles: First, the boundary con- 
ditions ought to make Eq. ^TE\i solvable; second, the 
time-dependent boundary conditions must coincide with 
stationary boundary conditions ^ once the laser field 



is off (A(i) = 0). The following boundary conditions 
indeed fulfill these constrains 



c,{t) [1/r, -a, + 0(r,)], J = 1,...,A^,(29) 



where Cj{t) are unknown functions of time, and aj are 
characteristics of ZRPs, as in the field-free case. There- 
fore, we have reduce the problem of obtaining the wave 
function ^' to the problem of calculating the functions 
Cj{t). These functions can be found once we impose the 
wave function given by Eq. (j26[) to the boundary condi- 
tions P5)) . Indeed, calculating the limits rj — > in Eq. 
(PS)) and using Eqs. (P^)) and (|B4|) . we achieve the system 
of integral equations for fj(t) {j = 1, . . . , N): 

N 

.,-V2i) =^Gfe(R,) 

k = l 



dxe" 



[W{t,x)f,{t~x)~f,{t)], (30) 



(ix)3/2 

where W{t,x) is given by Eq. (|B5p . The functions 
Cj{t) and fj{t) are connected by the equality: Cj{t) = 
~i^/2TT exp{—iet) fj (t) . 

We summarize results obtained until now in the form 
of a definition: Photoionization of a one electron system 
bounded in the field of N ZRPs, which are characterized 
by aj, is defined by the wave function ^'(r, t) 



*(r,<) = exp(^ 



N 

-*6i)VG,(r), 



(31) 



where Gj (r) is given by Eq. (jST)) , and the functions fj (t) 
and the parameter e are determined by solving non- linear 
eigenfunction-eigenvalue problem (|30p . 

However, solving non-linear eigenfunction-eigenvalue 
problems is generally a challenging task (for review see, 
e.g., Ref. [9^1). In the case of monochromatic radiation, 
there exists the Floquet twin-transformation (see, e.g., 
Refs. [20l. [97}). which tremendously simplifies numerical 
solution of problem ([30| . 

Note that all the consideration has been done only in 
the velocity gauge. Therefore, it is vital to check whether 
such a treatment is gauge invariant. It is well known that 
the wave function 'I'(r,t) is transformed as 

^'(r, t) ^ *(r, t) - cxp[^x(r, t)]^{r, t) (32) 

under some gauge transformation. Substituting Eq. (j29[) 
into Eq. ()32|) . we obtain the boundary conditions for 
^'(r,t): 



l-., + .£i:^^M + G(.,) 



(33) 



where Cj[t) = exp[zx(Rj, t)]cj (t). Having applied our 
renormalization procedure to this new wave function, we 
obtain 



N 



*(r, t) = exp[-zei + ix(i-, i)] 51 (r) 



(34) 



8 



where Gj{t) are given by Eq. (|27|) after performing the 
substitution fj{t) fj{t) = exp[—ix{'R-j,t)]fj{t); the 
functions fj{t) satisfy Eq. ([30]) where fj{t) — > fj{t) 
and Gj{t) Gj{t); as in the velocity gauge, Cj{t) = 
— i-\/27rexp(— iei)/j(t). Roughly speaking, the gauge 
transformation factor is "absorbed" by fj{t) such that 
all the equations are left form invariant. Hence, we have 
explicitly demonstrated that the developed technique is 
indeed gauge independent; as a result, it does not "suf- 
fer" from the "curse of the displaced atom," which is an 
artifact associated with the choice of gauge in the strong 
field approximation [95| . 

We present the special case of Eqs. ([50]) and PT|) for 
a single ZRP [N = 1) located at the origin (Ri — 0): 
The integral representation of the wave function ^'i(r, t) 
reads 

vl/i(r,i) = (27r)3/2 /* e-"*'i^v-G(r,i|0,i')/(i')rfi',(35) 



where the unknown function, /(t), and parameter, e, 
obey the integral equation 



i^K. + i\/2€j f{t) 



,3/2 



[W{t,x)f{t-x)- f{t)]. (36) 



This solution resembles the previous results [T^, [T^ in 
the case of a monochromatic laser field. 

A vital step in our derivation was the Floquet ansiitze 
(fT9|) . Note that in the case of ionization by a periodic 
laser field, this ansatze is the result of the Floquet the- 
orem, where e is called the quasi energy and its imag- 
inary part is proportional to the total ionization rates, 
and $(r, t) being a time periodic function. In the case of 
a periodic field, the Floquet theorem as well as quasi en- 
ergy formalism is extensively used in str ong field physics, 
for reviews see, e.g., Refs. [H HI HI Eo^- Therefore, 
the Floquet ansatze can be regarded as a trick to extend 
methods developed for periodic laser fields to the general 
case of arbitrary laser pulses. 



B. Stimulated Bremsstrahlung 

Once we have worked out photoionization, the problem 
of stimulated bremsstrahlung is trouble-free. In some 
sense, bremsstrahlung is even easer because there is no 
need to employ the Floquet ansatze [Eq. (fT9)) ]. 

We want to find a solution of Eq. ([T8|) with the 
initial condition being an outgoing-wave wave function 
(fO)l . iV'b)- Let \^i){t)) denote such a solution, and JJ 
represents the total propagator of the system at hand, 
= U{t, —oo) IV'b). From the initial condition, we 
conclude that 



lim UvG{t,to)\i^b) 

tn — ' — oo 



= (2^)'^/2exp|-^ J\p + A{r)fdr^ |p) . (37) 

Substituting Eq. (|37p into the Lippmann-Schwinger 
equation: 

Uit,to) = UvGit,to)~i f UvGit,t')VUit',to)dt', 



the following integral representation of the total wave 
function is obtained 

^'6(r,t) =exp|-^ J [p + A{T)f dr + ip ■ rj 
H'^^Y'^Y. / KvG{vA'^j.t')q,{t')dt', (38) 



where qj{t) = i^j*fc(Rj, t)/(27r)''/^. In the three- 
dimensional case [d — 3), the renormalized functions 
qj{t) obey the following system of non- homogeneous in- 
tegral equations 

iV^najqjit) = cxp J [P + dr + ip ■ Rj| 

-(2^)3/2^ / KvG(£^3,t\KkX)qk{t')dt' 

k = l J -oo 

[W{t,x)q,{t-x)-q,{t)]. (39) 



r°° dx 



APPENDIX A: THE RENORMALIZATION 
TECHNIQUE FOR BOUND STATES OF 
NONZERO ANGULAR MOMENTUM 

In this appendix we restrict ourself to the case of 
there-dimensions {d = 3) and a single ZRP centered at 
the origin (Rj = 0). Then, the heuristic potential ([9]), 
(ri| V |r2) = — /i(5(ri — r2)(5(ri), can be rewritten in the 
following form 



(k\V\r) = -AiV4^yoo(f^k)<5(r 



(Al) 



where ^(^(^k) is the spherical harmonic, fit denote the 
pair of the spherical angular coordinates of the vector 
k. Generalizing Eq. (jAip . we define a new heuristic 
potential 



(A2) 



(k| Vi„, |r) = -tiVATrk%rnink)Sir 



Now, we solve the stationary Schrodinger equation with 
this potential. 



py2 + Vi„A \^)^E\ij). 



The bound state solution reads 

(k| Vi„, 1^) k^YU^ 
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where (3 = 2-\/47rA*'0(O) is similar to Eq. (fTTI) . Fourier 
transforming the wave function in the momentum repre- 
sentation, e.g., by employing the partial wave expansion, 

exp(^k• r) = 4^ Ji'ikr)Y;^,{nu)Yv„Anr), 

V , m' 

we obtain 

V'(r) = B '— Yira{^r), (A4) 



One can readily notice that the expression under the sec- 
ond integral in Eq. (|Bip is regular at a; = 0. Using the 
following integral representation of the Macdonald func- 
tion [li 



poo 




iz 




/ exp 




2" 


X 


Jo 


A ^ ) 





dx. 



such that Im (z) > and Im (a^z) > 0, we obtain 



where i? is a constant proportional to (3. The wave func- 
tion (jA4p also represents the wave function of a particle 
with an arbitrary angular momentum I outside the spher- 
ical potential well. Evidently, Eq. (|A4[) coincides with 
Eq. ^ for an s-state, I = 0. 



APPENDIX B: THE ASYMPTOTIC BEHAVIOR 
OF THE FUNCTION G^(r) FOR SMALL 

In this appendix we calculate the asymptotics of the 

0. 



function Gj{r) [Eq. ^] when r ^ Rj, i.e. 
Since an arbitrary propagator K (r, t\r' , t') obeys the fol- 
lowing normalization condition 

YiinK{Y,t\v' ,t') = (5(r-r'), 

the function Gj{r) for any propagator K may have a 
singularity at r-,- — > 0. 

Let us separate this singularity. Having changed the 
variable of integration in Eq. (I??)) to a; = t — t', we can 
represent Gj(r) as 



dx e'^^ 



°° dx 



-^^^ A2(r)dr|/,(t-a;)-/,(i)cxp|^^^ . (] 



dx / ir^ ^ 
— -r- exp — - + lex 



= 2(^^j (e— /^v^r,) , (B2) 

where the following convention has been adopted 

z'^ = |z|'^exp(icargz), — tt < argz < tt. (B3) 



Substituting Eq. (|B2|) into Eq. jBT]), wc find the sought 
asymptotic in the three-dimensional case 



G,(r) 



[W{t,x)f,{t-x)~f,{t)] 



Jo (ix)^/^ 
+0{r,), (d = 3) 

where the function W(t, x) is given by 

W{t,x) = exp|^(^^ A(r)dr 



(B4) 



A^(T)dT , (B5) 
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